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Abstract 

We characterize a novel orthorhombic phase (7) of NaAlH4, discovered using first-principles 
molecular dynamics, and discuss its relevance to the dehydrogenation mechanism. This phase 
is close in energy to the known low-temperature structure and becomes the stabler phase above 
320 K, thanks to a larger vibrational entropy associated with AIH4 rotational modes. The struc- 
tural similarity of 7-NaAlH4 to o-NasAlHg suggests it acts as a key intermediate during hydrogen 
release. Findings are consistent with recent experiments recording an unknown phase during de- 
hydrogenation. 
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Complex light-metal hydrides represent a materials solution to the hydrogen-storage 
problem [if because of their potential for acceptably high gravimetric capacities, as 
well as an ability to be engineered for favorable hydrogen release kinetics. Of these, 
sodium alanate (NaAlH 4 ) has attracted particular interest by combining relatively high 
theoretical hydrogen release with ready reversibility through the addition of transition 
metal dopants jj B]. As such, the material has been widely studied as a template 
for developing novel hydro gen storage solutions. Several recent first-principles stud- 
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have addressed the structure and 



thermodynamics of the sodium alanate system. Nevertheless, important questions regarding 
the precise mechanism of dehydrogenation and the nature of the key phase transition from 
NaAlH 4 to NasAlHg remain unanswered. 

In this Letter, we present static, linear-response, and molecular dynamics (MD) calcu- 
lations based on density-functional theory (DFT) that support the existence of a hitherto 
unknown orthorhombic phase of NaAlH 4 (7), found to be thermodynamically favored over 
the low-temperature a phase at relevant operating temperatures for dehydrogenation. This 
phase is entropically stabilized by activation of AIH4 rotation and readily nucleates in MD 
simulations of (001) surface slabs. Our evidence is presented alongside recent experiments to 
highlight the likely significance of this phase in stabilizing and mediating the decomposition 
reaction. 

Our calculations are based on DFT in the plane- wave pseudopotential formalism 21] 
using the PBE exchange-correlation functional 22|. Ultrasoft pseudopotentials 23j for Na 
(2s 2 2j? 6 3s a5 ), H, and Al (3s 2 3p 2 ) were taken from the Quantum-Espresso distribution 21]. 
A second Troullier-Martins pseudopotential [24] for Na, used only for the MD runs, was 
generated in a 3s°' 5 3p 05 valence configuration, with nonlinear core corrections [3] added. 
For the MD simulations, cutoffs of 25 and 200 Ry were used for the wavefunctions and 
charge density, respectively, representing convergence of forces to within 0.02 eV/A. These 
were raised to 30 and 300 Ry for the total-energy and linear-response calculations. The 
a — > 7 transition was studied using Car-Parrinello MD simulations (26] of bulk supercells of 
NaAlH 4 with 96 atoms (2x2x1 a unit cells) and of slab supercells with 192 atoms (8 layers; 
2x2x2 a unit cells). Slabs were cleaved along the (001) plane — verified experimentally 
and theoretically as the stablest facet [l^] — with 14 A of vacuum inserted between images. 
Except where indicated, slab simulations were performed in the NVT ensemble, whereas 
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bulk simulations were performed in the NPT ensemble at zero pressure using the Parrinello- 
Rahman extended Lagrangian formalism [27]. The theoretical a lattice parameter was used 
for the initial cell state. Bulk NaAlH 4 was sampled at 50 K intervals from 300-650 K; surface 
slabs at 25 K intervals from 225-425 K. Temperatures were maintained using Nose-Hoover 
chains 28], with electronic fictitious mass \x = 500 au and At = 6 au. Following 5 ps of 



thermalization, each bulk simulation ran for 25 ps and each surface simulation for 15 ps. 

The new phase, designated throughout this Letter as 7 to distinguish it from the high- 
pressure j3 phase [29|, manifests itself in our MD simulations of the (001) surface slab of 
NaAlHzi. Specifically, we find that the (001) slab exhibits a spontaneous transition from the 
a to the 7 phase for 250 < T < 350 K, within 15 ps into the production run. This transition 
is clearly visible in Fig. [T^i, which compares final configurations of the surface slab at select 
simulation temperatures. (Melting occurred above 350 K; tendency to melt is enhanced by 
our finite slab thickness). To rule out finite-size or supercell restrictions, we ran additional 
surface simulations at 300 K, one with a larger 270-atom supercell and one in the NPT 
ensemble. Both evidenced the same transition to the 7 phase. An NPT simulation of bulk 
7-NaAlH4 at 300 K was also performed, with the new phase remaining stable throughout 
the run. 

Figure [lb illustrates the basic structural differences between the a and 7 phases, as 
observed in the slab simulations. Three qualitative features of the a —> 7 transition emerge 
in the MD runs: a rotational disordering of the AIH4 groups; an expansion of the crystal 
along c; and a shear of successive crystal planes parallel to the surface, generating the lattice 
symmetry of the 7 phase. 

In contrast to the surface slab, no structural reorganization was observed in NPT simu- 
lations of bulk NaAlH4 (until melting at T > 550 K). Given its facile expression in the (001) 
slab, we conclude that any observation of the transition to 7 in the bulk must be kinetically 
hindered, requiring much longer simulations for observation. This hindrance relates to the 
volume expansion of the lattice during the transition, as discussed further below. 

In Fig. [21 we present the structure of 7-NaAlH4, taken from representative dynamics 
timesteps and relaxed with respect to all internal parameters and lattice vectors. The unit 
cell is best described as a 12-atom base-centered orthorhombic (bco) structure with space 
group Pmmn and a = 6.500 A, b = 7.105 A, and c = 7.072 A (see Table [I] for details; note 
the near-degeneracy of b and c). We compared the energy of this structure, calculated using 
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FIG. 1: (Color online) (a) Equilibrated structure of an 8-layer (001) surface slab of NaAlH4 taken 
from MD runs at 225, 250, and 275 K (topmost four layers shown). Na/Al/H atoms are shown 
in blue/green/white. Slabs are viewed along the b axis of the a phase, with c oriented upwards. 
Surface shear at 250 K signals the onset of 7, complete by 275 K. (b) Schematic illustration of 
structural differences before (left) and after (right) the transition (identical viewpoints, Na atoms 
omitted). Solid lines outline the a and 7 unit cells; dotted lines connect original boundary atoms 
of the a cell. 

216 k-points in the full Brillouin Zone, to that of the low-temperature a phase, calculated 
using a 324-k point mesh of similar density. At K, the a phase is lower in energy by only 
9 meV/atom, pointing to a close competition for stability between these two phases. The 
accuracy of the PBE energetics was also checked against quantum-chemistry calculations on 
isolated 72-atom clusters of a- and 7-NaAlH 4 using the B3LYP functional and a 6-311+G 
basis; the difference between B3LYP and PBE in the relative stability of the two clusters 
(likely bracketing the exact result) was 2 meV/atom. 

A comparison of Figs.[T]and[2]reveals differences between the 7 phase at finite temperature 
and its zero-temperature, fully relaxed equivalent. The zero-temperature structure features 
alignment of H (but not Al) in successive A1H 4 tetrahedra along c. At finite temperatures, 
7-NaAlH 4 instead exhibits an average alignment of Al along c, a consequence of randomized 
AIH4 orientations and increased configurational freedom. 

To establish finite-temperature properties, we have calculated the phonons [30( and 
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FIG. 2: (Color online) Structure of orthorhombic 7-NaAlH4, shown (a) in the standard view; (b) 
along a; and (c) along c. The color scheme follows Fig. [TJ 

TABLE I: Internal atomic positions for 7-NaAlH 4 (Pmmn; a = 6.500 A, b = 7.105 A, c = 7.072 A). 
Units are fractional coordinates (a, f3, 7) of the primitive bco lattice vectors. 



Atom type 


a 


P 


7 


Na 


±0.141 


T0.141 


±0.250 


Al 


±0.349 


=F0.349 


T0.250 


H 


±0.198 


=F0.198 


T0.067 


H 


±0.198 


=F0.198 


T0.433 


H 


±0.313 


±0.317 


T0.250 


H 


±0.317 


±0.313 


±0.250 



phonon density of states (PDOS) for the a and 7 phases using density-functional pertur- 
bation theory 3l|. Figure [3] compares the PDOS for a-NaA!H 4 (in close agreement with 



Ref. Il6l) with 7-NaAlH 4 . Bands below 180 cm -1 are Na/Al translational lattice modes; 180- 
580 cm -1 , AIH4 rotational modes; 750-900 cm -1 , bending modes; and above 1600 cm -1 , 
stretching modes. The most essential difference in the phonon spectra is that the AIH4 ro- 
tational modes (350-580 cm -1 in the a phase) are softened substantially (200-420 cm -1 ) in 
the 7 phase. This translates to an increase in vibrational entropy at temperatures for which 
these modes are active. We estimate the temperature-dependent free energy F(T) of the 
phases by adding the Bose-Einstein vibrational entropy contribution onto the ground-state 
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FIG. 3: Phonon density of states (grey shaded region) and Raman spectrum (solid black line, with 
thermal broadening) for 7- (top) and a-NaAlH4 (bottom). Raman peaks below 1000 cm -1 are 
amplified tenfold for easier visibility. The arrow marks the the a — » 7 transition temperature. 
Inset shows the temperature dependence of the free-energy difference AF = — F a . 



total energy E : 



F(T) = E + J2 ll^U + k BT\n 



1:1 



exp 



q..y 



where sums run over band indices j and wavevectors q in the first Brillouin Zone. Neglecting 
further contributions to F(T), we find good agreement between the transition temperature 
predicted by the free-energy analysis and that found in the dynamics simulations (250 K) 
when using identical simulation parameters. Upon increasing the plane-wave cutoff and 
including the Na semicore states to more accurately account for the zero-point motion as 
reported in Ref. [la, we obtain the results in Fig. [3j indicating a thermodynamic preference 
for the 7 phase for T > 320 K (see inset). Notably, this establishes the 7 phase as the 
stablest variant within the experimental temperature range for dehydrogenation. 



In a recent in-situ Raman study 



32] , Yukawa et al. found evidence of an unknown inter- 



mediate phase of NaAlH-4 in first-stage dehydrogenation, emerging in the uncatalyzed sample 
just before melting at 458 K. The authors suggest that, like our proposed 7-NaAlH 4 , this 
phase should be stable upon cooling to near room temperature and contain intact AIH4 units 
with little constraint from the surrounding crystal lattice. Accordingly, we have calculated 
the Raman spectrum for 7-NaAlH 4 , shown in Fig. [3] 3y]. A comparison between our results 
and the spectrum at 458 K in Ref. [32| shows a remarkable agreement between the two. Fore- 
most, we note the new signature Al-H stretching peak near 1800 cm -1 in both experiment 
and simulation. The simulations also record a second resonance at 1750 cm -1 ; although not 
explicitly visible in Ref. |32j, asymmetry in their 1800 cm 1 signature does suggest the pres- 
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ence of an additional lower-frequency peak. The difficulty in experimentally resolving these 
two peaks can be attributed to the thermally averaged nature of the finite-temperature 7 
structure with respect to its zero-temperature counterpart, as well as additional soft lattice 
modes activated during the transition (similar broadening is observable in the pair corre- 
lation data of Fig. H]). In the experiment, the 400-500 cm -1 peaks disappear, and distinct 
peaks at 770 and 817 cm -1 merge into a single broad peak. Both features are also visible 
in Fig. [31 Note that in reality, the 400-500 cm -1 rotational modes in a-NaAlH4 do not 
disappear in 7 but rather shift to a broad peak at ~300 cm -1 ; however, the corresponding 
Raman-active peaks are too weak to demonstrate a resolvable signal. Collectively, these 
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similarities present a strong case for the interpretation of the unknown phase in Ref. 
7-NaAlH 4 . 

A closer examination of the MD results elucidates why the a — » 7 transition is not ob- 
served in our bulk NPT simulations. In the surface slab simulations, 7-NaAlH 4 nucleated 
only upon activation of the rotational modes of surface AIH4 units (center panel, Fig. UK). 
This agrees with Fig. [31 since softening of AIH4 rotational modes is the predominant contrib- 
utor to the thermodynamic stability of the 7 phase. However, in bulk a-NaAlH4 the AIH4 
units cannot rotate: the barrier to rotational mobility is larger than for the relatively uncon- 
strained surface groups, and rotation must be coupled with large volume fluctuations. This 
makes ready nucleation of the 7 phase in bulk NaAlH 4 difficult. Correspondingly, rotation 
of interior AIH4 units in the (001) surface slab is strongly coupled to the volume expansion 
perpendicular to the surface. The calculated surface energy of 7-NaAlH 4 (001) is also much 
smaller than a-NaAlH^OOl) (3.1 versus 8.2 meV/A 2 ), favoring the transition in the slab. 
The results show that A1H 4 rotation percolates once activated at the surface, as the slab 
thickness increases to accommodate. This lowers the barrier for afe-planar shear, finalizing 
the transition. We therefore assume that nucleation of 7-NaAlH 4 should occur more readily 
at an exposed (001) surface, at grain boundaries in a polycrystalline sample, or in low-stress 
regions that can accommodate the local volume fluctuations necessary for AIH4 rotation. 

We investigate the transition timescale in the MD simulations by introducing an order 
parameter A as a reaction coordinate, defined as A = 1 — /edge //comer- Here / e d ge and /comer 
are the fractional occupancies of the edge-center and corner sublattices in the projection 
of the a unit-cell geometry onto the surface plane (see Fig. [lb) . The rapid evolution from 
A = (i.e., pure a) to A = 1 (i.e., pure 7) evidences the strong thermodynamic driving force 
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FIG. 4: Comparison of the Al-Al pair correlation function <?ai-A1 f° r the (001) slab before (dashed 
line) and after (solid line) the a — > 7 transition at 275 K. Solid vertical lines give peak locations 
for a lattice-matched ideal crystal of a-NaaAlHg (see text for definition), with relative peak heights 
shown for comparison. Inset shows the evolution of the order parameter A during the transition. 

for the transition (see inset, Fig. H]). 

Geometric similarities in the aluminum substructures of 7-NaAlH 4 and a-Na3AlH 6 sug- 
gest a phenomenological connection between the a — > 7 transition and the NaAllHLi — > 
NasAlHg dehydrogenation reaction. These are quantified in Fig. HI which compares the 
Al-Al pair correlation function for the (001) slab during the a — > 7 transition with one for 
Na3AlH~6 that is constrained to be lattice matched to the 7 cell. Specifically, a-NasAlHg 
peaks were obtained by isotropically shrinking the cell volume (-10%) to match the mean 
lattice parameter in the ab plane to the 7-NaAlH 4 cell boundaries. The striking agreement 
between the locations and relative heights of the two sets of peaks strongly hints at the 
role of 7-NaAlH~4 as a key intermediate during dehydrogenation. If so, the role of exposed 
surfaces in nucleating 7 may help to account for the enhancement of dehydrogenation ki- 



netics due to ba 
nanoparticles 



1 milling without a catalyst [33], as well as improved kinetics of NaAlH4 
34|. We did not examine the impact of transition metal additives; however, 
the catalyst is likely to further aid the conversion of 7-NaAlH 4 to a-Na3AlH 6 . A detailed 
investigation of this possibility is currently underway. 

In conclusion, we have discovered and characterized a new phase of NaAlrFi, labeled here 
as 7. Our description is fully consistent with an unknown intermediate found in a recent 
in-situ Raman study. Activation of the A1H 4 rotational modes drives the transition to 7, 
which occurs readily in (001) surface slab simulations for T > 320 K. The thermodynamic 
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stability of 7-NaAlH4 beyond this temperature has been further confirmed by vibrational 
entropy calculations using density-functional perturbation theory. Remarkable structural 
similarities between 7-NaAlH 4 and a-Na3AlH 6 suggest that the 7 phase should play a key 
role in mediating dehydrogenation. Notably, no Al-H bonds were broken in our simulations; 
this agrees with Ref. |35| and underscores the need for an additional defect-driven mechanism 
to fully account for dehydrogenation. Accordingly, further investigation of low-energy defects 
in 7-NaAlH4 that could facilitate Al-H cleavage is recommended. 
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Calculations were performed using Quantum-Espresso [2jJ on facilities provided through 
NSF Grant DMR-0414849. The authors thank Gerbrand Ceder for helpful discussions, and 
Elise Li for the quantum-chemistry results. 
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